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Between 1995 and 2009, electron temperature (T e ) measurements of more than 15000 plasmas produced in 
the Joint European Torus (JET) have been carefully reviewed using the two main diagnostics available over 
this time period: Michelson interferometer and Thomson scattering systems. Long term stability of JET T e 
is experimentaly observed by defining the ECE TS ratio as the ratio of central T e measured by MICH and 
LIDAR. 

This paper, based on a careful review of T e measurement from 15 years of JET plasmas, concludes that 
JET T e exhibits a 15-20% effective uncertainty mostly made of large-scale temporal drifts, and an overall 
uncertainty of 16-22%. 

Variations of 18 plasma parameters are checked in another data set, made of a reference data set made of 
ohmic pulses as similar as possible between 1998 and 2009. Time drifts of ECE TS ratios appear to be mostly 
disconnected from the variations observed on these 18 plasma parameters, except for the very low amplitude 
variations of Bo which are well correlated with off-plasma variations of a 8-channel integrator module used 
for measuring many magnetic signals from JET. 

From mid-2002 to 2009, temporal drifts of ECE TS ratios are regarded as calibration drifts possibly caused 
by unexpected sensitivity to unknown parameters; the external temperature on JET site might be the best 
parameter suspected so far. 

Off-plasma monitoring of MICH made of calibration performed in the laboratory are reported and do not 
appear to be clearly correlated with drifts of ECE TS ratio and variations of KC1D integrators. Comparison of 
estimations of plasma thermal energy for purely Ohmic and NBI-only plasmas does not provide any definite 
information on the accuracy of MICH or LIDAR measurements. 

Solutions aiming at tracking down these unexpected uncertainties of JET T e are detailed and can be 
performed during next JET campaigns (C28+, after October 2011), for instance with highly-reproducible 
reference pulses and off-plasma monitoring of the diagnostics. 

Whatever causes these T e drifts, this experimental issue is regarded as crucial for JET data quality. 



I. INTRODUCTION 

In this paper 15 years of electron temperature (T e ) 
measurements made by multiple diagnostics on plasmas 
performed in the Joint European Torus (JET) have been 
analysed and carefully reviewed. 

JET offers the unique possibility to assess T e stabil- 
ity and, then, its effective uncertainty, by checking the 
very long-term stability of electron temperature measure- 
ments provided by multiple diagnostics on a large sized 
tokamak. Lessons learn from this data review are then 
expected to provide solid feedback for JET operation and 
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benefit future ITER operation. 

In II A, diagnostics measuring the electron tempera- 
ture profile on JET are described: Michelson interfer- 
ometer (MICH), core Thomson scattering (LIDAR), ECE 
heterodyne radiometer (HRAD) and the high resolution 
Thomson scattering system (HRTS). Some details about 
the absolute calibration of the ECE diagnostics (MICH, 
HRAD) are given in (II B). More than fifteen years of 
JET data are analysed in this paper, with different time 
and space resolution and non-negligible data scattering 
(II C). The way T e data from different diagnostics are 
compared is stated in (II D) for the magnetic mapping 
of the diagnostics' line of sight. JET T e observables are 
defined in (HE). 

In III, experimental observations of the long term sta- 
bility of JET T e between 1995 and 2009 are done by 
using the ECE TS ratio (III A). Two data sets are used 
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in this paper (IIIB): the large data set, from 1995 to 
2009, where MICH and LIDAR data are available and 
with a broad selection of plasma parameters (IIIC); and 
the reference data set made of ohmic pulses as similar as 
possible between 1998 and 2009, including T e data from 
the four diagnostics (HID). Reduced set limited to the 
latest JET campaigns is used as well, ranging from cam- 
paign C20 (April 2008) to C27 (November 2009) i.e. JET 
pulse numbers (JPN) ranging from 72639 to 79853. 

In IV, the evolution patterns and/or instrumental 
drifts of ECE TS ratio are linked with small variations 
of the magnetic field, not caused by or linked to JET 
plasmas, namely off-plasma. 

In V, the focus is set on MICH as off-plasma monitoring 
has been carried out in the laboratory since 1984 - this is 
used to assess the stability of the diagnostic. Data from 
1998 to 2010 are presented and compared with the time 
drifts of the ECE TS ratio. 

In VI, the plasma thermal energy (W t h) is estimated 
from a kinetic expression, that is using T e profiles mea- 
sured by MICH and LIDAR, then compared to W t h from 
diamagnetic measurements. 

In VII, the profile mismatches of high resolution mea- 
surements from ECE (HRAD) and TS (HRTS) systems 
are documented. Links are made with ECE TS ratio. 

The findings of the study are then discussed in (VIII), 
where the overall uncertainty of JET T e measurements 
between 1995 and 2009 is assessed. Comments and rec- 
ommendations are then provided and solutions to track 
down the unexpected variations of the ECE TS ratio are 
detailed. 



II. DIAGNOSTICS & DATA 
A. Diagnostics 

In the JET tokamak, radial profiles of electron tem- 
perature are measured by diagnostics based on electon 
cyclotron emission (ECE, Refs 7 ' 14 ' 25 ): Michelson inter- 
ferometer (MICH), heterodyne radiometer (HRAD); and 
based on Thomson scattering (TS, Refs 14 ' 25 ): LIDAR 
and a high-resolution TS system (HRTS). MICH and LI- 
DAR are both absolutely calibrated while HRAD is cross- 
calibrated. Even if HRTS could be independently cali- 
brated, HRTS data used in this article are actually cross- 
calibrated. Main characteristics of these diagnostics are 
summarized in Table I. 

The Michelson interferometer is originally described 
in Ref 6 in its original 1983 set-up; a description of 
the current (2012) set-up can be found in Ref 23 . The 
absolute calibration of MICH used in this study has 
been performed in 1996 using the usual cold/hot sources 
technique 1 ' 4 based on Ray leigh- Jeans' law. Such a cal- 
ibration requires the positioning of sources in front of 
the diagnostic in-vessel antenna: the interferometer can 
then be calibrated against these sources, characterised 
by their surface temperature, emissivity and frequency 



domain. MICH is calibrated for the 50 GHz to 350 GHz 
frequency domain, with a 10% uncertainty. 

The HRAD diagnostic is a heterodyne radiometer in- 
stalled in 1987 21 in a 8-channel configuration cover- 
ing the 73-79 GHz range. Constantly upgraded (see 
Refs 2 ' 3 ' 17 and references herein), 96 channels are available 
since 2008, covering 69-207 GHz range. HRAD is cross- 
calibrated against MICH for each JET pulse by equalising 
signals from each radiometer channel to the temperature 
measured by MICH for the same ECE frequency 27 . A 
fully calibrated radiometer is technically not out of reach 
but absolute calibrations should be performed after each 
modification or upgrade. Furthermore high-frequency 
microwave components (MOO GHz) are quite prone to 
failure, strongly reducing the effectiveness of such a cali- 
brated radiometer. 

The LIDAR 12 is a Thomson Scattering diagnostic 
based on a time-of-fiight technique, first installed in 1987, 
where radial position is determined by the time it takes 
the scattered light to reach the spectrometer. In the 
spectrometer the incoming light signal is divided into 
several spectral intervals (370 nm to 850 nm) by inter- 
ference filters and then focused into fast photosensitive 
detectors. Calibration of the system is performed in sev- 
eral steps. First the spectral response of each individ- 
ual channel is measured by illuminating the spectrom- 
eter input slit with a monochromatic light source with 
adjustable wavelength. The whole collection system is 
illuminated by a calibrated white light source placed in 
front of the JET vacuum windows and the relative sensi- 
tivity of each channel is recorded. Spectral transmission 
of the windows is measured separately and applied as 
a correction to the calibration factors. Window trans- 
mission is chromatic (depends on wavelength) and the 
associated correction does cause a change in measured 
T e but not more than 5%. Window transmission was 
measured several times during 2005-2011 and appeared 
to not change significantly. The spectral calibration is 
relatively easy to perform and tends to show the same 
results every time, and as has been stable for many years 
(> 15). Systematic error in T e measurements caused 
by these calibration uncertainties is very unlikely to ex- 
ceed 5%. Position calibration is performed by observing 
the light scattered from the laser beam dump with well 
known location at the JET inner wall. 

The HRTS diagnostic 20 measures the electron temper- 
ature (T e ) and density (n e ) along a chord that runs near 
the plasma mid-plane, from the plasma core to the outer 
pedestal region. This profile is measured by analysing 
the light that is Thomson scattered from a 3 J, 20 Hz, 
Q-switched NdiYAG system operating at 1064 nm. The 
scattered light is coupled into a linear array of more than 
130 optical fibers which feed into a bank of 21 polychro- 
mators. Each polychromator analyzes three separate sig- 
nals, via delay-line multiplexing, allowing a maximum of 
63 spatial points in the HRTS profile. The spectral re- 
sponse and absolute transmission level for each spatial 
point is independently calibrated. The absolute position 
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FIG. 1. Line-of-sights of MICH, HRAD, LIDAR and HRTS 
for typical ohmic plasma (JPN 79853, t=5.0s). Flux surfaces 
are computed by EFIT. 

of the profile is performed by back-illuminating the opti- 
cal fiber array onto an in-vessel ruler. Due to (corrected 
for data acquired after 2010) problems with the spectral 
calibration process, the HRTS T e data from JET cam- 
paigns for C23-C27 (JPN =74391-79853) relies in part 
on cross-calibration with the HRAD diagnostic 10 . This 
cross-calibration is performed by comparing the T e pro- 
file with the HRAD diagnostic for a handful of carefully 
selected dedicated calibration pulses. The resultant cali- 
bration correction factors are then validated by compar- 
ing the corrected profiles over many 1000's of JET pulses. 



Diagnostic 


measurement 
time 


measurement 
frequency 


radial 
definition 


MICH 


17 ms 


30 Hz 


-10 cm 


HRAD 


200 us 


5 kHz 


2-5 cm 


LIDAR 


5 ns 


4 Hz 


-12 cm 


HRTS 


20 ns 


20 Hz 


-1 cm 



TABLE I. Temporal & radial characteristics of JET diagnos- 
tics measuring T e 

The diagnostics' lines-of-sight (LOS) are plotted in 
Fig. 1. For vessel configurations used between 1995 and 
2009, MICH and LIDAR LOS are usually very close to the 
plasma midplane by a few cm. HRAD and HRTS LOS are 
situated at lower z-position and usually miss the plasma 
centre. Mapped on the plasma midplane, HRAD T e pro- 
files usually exhibit a gap for \p\ < 0.2 and HRTS ones 
are restricted to the low field side (p > 0.2 — 0.3), where 
p = an d 4> the normalised poloidal flux. 

B. About ECE calibrations 

MICH's absolute calibration is used by HRAD and 
HRTS (for C20-C27); these diagnostics provide the JET 
T e radial profiles with high time and radial resolutions 
that are usually shown in publications. 

Absolute calibration of MICH follows, for each fre- 
quency F: 

T keV (F)=I meas (F)/a(F) 



with imeas(^) the uncalibrated intensity measured by 
MICH and a(F) the F-dependent calibration factors. 

A radial calibration of the ECE system requires a 
magnetic reconstruction code (EFIT 15 ' 19 ) that links the 
frequency domain of the observed ECE to the corre- 
sponding radial domain. The usual assumption of a 
cold ECE resonance is made, which means that uj n = 
nq e B§jm e i.e. all broadening effects, including relativis- 
ts ones, are neglected; see Ref 5 for a detailed analysis. 
This assumption results in a slight (a few cm) artificial 
shift of the ECE T e profiles in the outwards direction 

(KECE,^cold < RECE,cold)- 

In order to update the (undocumented) 1996 ECE cali- 
bration in use for data measured between 1995 and 2010, 
a calibration campaign has been performed in 2007 26 and 
repeated for confirmation in 2010 22 , addressing most of 
the issues raised in the conclusion of the 2007 campaign. 

Both calibration campaigns indicate a small decrease 
of MICH's sensitivity in the relevant ECE range (90 GHz 
to 220 GHz), causing (if applied) a 15 — 20% increase 
of JET T e . Through cross-calibrations procedures (see 
(II A)), HRAD and HRTS measurement would be affected 
similarly. 

Following this revaluation, the periods of validity of 
the 1996, 2007 and 2010 ECE calibrations have been as- 
sessed. Data from MICH, HRAD and HRTS presented in 
this study still use the 1996 calibration. Applying the 
2010 ECE calibration leads to a noticeable increase of T e 
affecting the ECE TS ratio but not its variations. 



C. Data 

Data analysis mostly focuses on MICH and LIDAR 
measurements, which might be highly scattered: ±5% 
for MICH and ±10-15% for LIDAR. Time averaging is 
then required. JET pulses have been cut into timeslices 
of Is length where the total field (Bo), plasma current 
(Ip), additional power (Padd), electron temperature and 
density (T e , n e ) are kept reasonably constant: < 1% vari- 
ations for Bo, < 5% for Ip, < 10% for central density, dia- 
magnetic energy and additional power. All the plasma 
measurements presented in this paper are extracted from 
the JET public database. 

In order to improve the readabiliy of the figures, tem- 
poral axes often represent timeslices instead of pulse 
dates or JPN. Timeslices selection is roughly linear with 
JPN as pulses usually have 5 to 10 timeslices with ohmic 
conditions, and not linear with pulse date. This choice is 
relevant as long term variations of plasma measurements 
(usually T e ) are then easily seen. 

All JET pulses between 35000 <= JPN <= 79853 
are included in the forthcoming analyses, where JPN 
35000 was produced in May 1995 and JPN 79853 is the 
last pulse of C27 in November 2009. Pulse selection 
(pulses/timeslices) is done on plasma parameters as de- 
tailed in (IIIB), providing of course that pulses were not 
plasma-free and that the ECE and TS diagnostics were 
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operating. 

Routines performing the data processing for all di- 
agnostics, including ECE and TS as well as all other 
JET data, are constantly improved, which usually im- 
plies some data reprocessing. JET data used and shown 
in this study have been taken as stored in JET public 
database during Summer 2011 and no specific reprocess- 
ing was performed. 



D. Magnetic mapping 

Temperature profiles in tokamaks are four dimensional 
quantities (three radial, one temporal). The assump- 
tion of toroidal symmetry reduces the radial dimensions 
to two. Comparing measurements from different diag- 
nostics, each having a different LOS, requires to map 
their geometrical LOS to the midplane plasma radius (or 
any other one) by using a magnetic equilibrium model; 
EFIT 15 ' 19 is used in this paper with a magnetic recon- 
struction only constrained by JET magnetic coils. 

Discussing the accuracy of JET ECE radial position of 
all JET data studied in this paper is definitely out the 
scope of this paper. The hypothesis of cold resonance 
brings a slight overestimation of ECE radial position - 
a few cm for usual JET plasma conditions. Mapping 
of LOS to the central midplane affects it as well; more 
generally, JET ECE radial position relies on magnetics 
measurements. 

Accuracy of JET magnetic measurements is 1.5% 1 . 
Following Ref 18 it leads to an accuracy of 1 cm to 1.5 cm 
for the magnetic geometric boundaries (for p = 1) as 
computed by EFIT. Magnetics measurements are not yet 
constraining EFIT enough when reconstructing the mag- 
netic equilibrium's radial position for central regions of 
the plasmas, radial accuracy at the plasma core is not 
expected 2 to be better than ~5-10cm. 

Radial localisation of ECE measurement depends on 
Bo, following Ar/r = AB/B for cold resonance hypoth- 
esis, r being the plasma major radius. This B -related 
uncertainty leads to Sr ~5 cm in the plasma centre. Tak- 
ing only into account this (quite large) source of radial 
uncertainty stresses the complexity of the problem. 

It is expected that these large uncertainties affecting 
the estimation of JET ECE radial position will not show 
clear trends of drifts with time, and that they are more 
likely to have the same effect as a low-amplitude white 
noise on T e observables (see (HE)). 

Measurement position of LIDAR and HRTS are cali- 
brated during JET shutdown and do not rely on EFIT on 
their LOS; they do though when mapping to the plasma 
midplane to be able to compare them to ECE T e profiles. 
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E. Comparing T e measurements 

1. T e observables 

Different quantities can be defined to reduce the bi- 
dimensional profiles to one single value: they are called 
observables and are designed to quantify the discrepan- 
cies between MICH and LIDAR. Some possibilities have 
been studied: 

• radial averaging around the plasma centre: 

T eMAX = median (T e (p,t)), \p\ < 0.2) 

• integrals of T e (p) over p where measurements were 
available for both diagnostics: 

/0.9 
T e(p)dp 
-0.4 

• T e values for fixed p values: T e0 .5, T e0 .8 and T e0 .g 
for, respectively, p = 0.5, 0.8 and 0.9. 

In order to reduce the uncertainties caused by a non- 
ideal midplane mapping (see (II D)), profiles are radially 
averaged over the plasma centre, following T eMAX . This 
method has been thoroughly tested and it is regarded as 
robust, satisfyingly dealing with high- or low-peaked pro- 
files as well as mapping issues around the plasma centre. 

Using T e _ .4^o.9 or Teo.5 instead of T eMAX would lead 
to very similar behaviour (variations and amplitudes). 
Mapping errors and measurement noise prevent T e0 .8 and 
T e0 .g to be give similar results. 

III. LONG-TERM OBSERVATIONS OF ECE TS RATIO 
A. Definition of ECE TS ratio 

Stability of JET electron temperature is based on ob- 
servations of T eMAX as measured by MICH and LIDAR, 
i.e. T eMAX ,MICH /TeMAx ,LIDAR; this ratio is called ECE 
TS ratio. 

The following paragraphs explain why HRAD and 
HRTS are not used to build an equivalent ECE TS ra- 
tio. 

HRAD data is calibrated for each pulse against MICH, 
defining cross-calibration factors T eMICH /T e 

HRAD ' 

any 

drift of MICH calibration could then be checked against 
these HRAD cross-calibration factors, provided the time 
stability of HRAD calibration is assessed independently. 
Despite our efforts, it has not been possible to use these 
HRAD cross-calibration factors to check MICH stability, 
as the former are obviously not stable enough to be used 
as references. Gaps in the cross-calibration factors can 
be explained by the waveguide switches that allow to se- 
lect O- or X- microwave mixers that actually have moving 
parts; large time drifts observed in these cross-calibration 



5 



factors (up to 96 for the latest set-up) do not all show 
similar features and each frequency channel appears to 
be almost independent from the other ones. 

HRAD should therefore be regarded in this paper as an 
upgrade of MICH, with higher time and radial resolution. 
Two points make MICH best suited than HRAD for the 
ECE TS ratio; first: MICH LOS is closer to the plasma 
centre than HRAD's; second: HRAD cross-calibration 
might add additional uncertainty when estimating T e . 

The main reasons to prefer LIDAR to HRTS to build 
the ECE TS ratio are, first, the late availability of HRTS 
(Spring 2008) and, then, the dependence of HRTS cali- 
bration to HRAD, then MICH, for all the data presented 
in this paper. 

For those reasons, ECE TS ratio is made of MICH and 
LIDAR data. 



B. About data sets 

Many plasma parameters can affect T e , the most obvi- 
ous probably being the amount of additional power cou- 
pled to the plasma. For non purely Ohmic pulses, a se- 
lection of pulses/timeslices has been made on Padd = 
Pnbi + PiCREb where Pnbi is the total additional power 
from NBI (Neutral Beam Injection) and Picrh is the 
total additional power from ICRH (Ion Cyclotron Radio- 
frequency Heating). Pulses with LHCD (Lower Hybrid 
Current Drive) are not included in the analysis because 
this heating system affects ECE signals, especially at the 
plasma edge. 

The following arbitrary separation on Padd = Pnbi + 
Picrh is then used: 

• Ohmic pulses with Padd =0 MW 

• "Low-power" pulses with Padd=0-10 MW 

• "High-power" pulses with P A dd=10-20MW 

Low-power pulses are mostly in L-mode while high- 
power ones are mostly in H-mode; no distinction is actu- 
ally made between both confinement modes in this paper. 

Two data sets are defined: 

• The "large data set" spans 15 years of data pro- 
duced in JET (1995-2009) with a broad selection 
of plasma parameters. This set gives an overview 
of the ECE TS ratio, i.e. how T e measurements 
from ECE (MICH) and TS (LIDAR) systems com- 
pare with time. 

• The "reference data set" spans 13 years of data 
produced in JET (1998-2009) with narrow selec- 
tions of plasma parameters (e.g. Bq = 2.41 T, Ip = 
1.96 MA); 18 parameters are checked and selected 
to be as constant as possible. These pulses are then 
regarded as "beacon" or "reference" pulses, allowing 
for more precise study of ECE TS variations with 
time. 



Pulse No: 35000-79853. ECE TS ratio. Ohmic data 




Values Values Values 



FIG. 2. ECE TS ratio for ohmic pulses, large data set, 1995- 
2009. 

(Top) Data in black, 50-timeslices average in red. Years are 
indicated 

(Bottom- left): distribution of ECE TS ratio data, with a half- 
width of 0.11 and Gaussian fit in red. 
(Bottom-centre): distribution for 50-timeslices average. 
(Bottom, right): distribution of residues, defined as ECE TS 
ratio minus averaged ones, with a half- width of 0.08. Gaussian 
fit in red 



C. Large data set: 1995-2009 
1. Description 

JET data has been analysed from 1995 to 2009 for 
35000 <= JPN <= 79853 where both MICH ad LIDAR 
were on-line. Data selection is limited on Bq and Ip with 
Bo limited to 2 to 3T and Ip to 1.5 to 4 MA; and after 
additional power separation: 

Ohmic: 15k pulses 135k timeslices 
Low-power: 6k pulses 23k timeslices 
High-power: 3k pulses 10k timeslices 



2. Observations: 1995-2009 

In Fig. 2, the ECE TS ratio is plotted for ohmic con- 
ditions. Raw data ranges from 0.5 to 1.25 and its dis- 
tribution is Gaussian-like (Fig. 2, bottom, left) with a 
mean around 0.95 and a half- width of 0.11. Averaged 
data are shown as well, over 50 timeslices that usually 
represents 5-10 pulses; distribution clearly departs from 
Gaussian (Fig. 2, bottom, centre) and residues look like 
white noise (null mean Gaussian-like distribution, Fig. 2, 
bottom, right) with a 0.08 half- width. Evolution patterns 
are obviously noticeable on raw and averaged data. 

Focusing on averaged data, ECE TS ratio shows very 
clear structures roughly centered on ~0.95 with ~20% 
temporal variations (peak-to-peak). 

Ohmic, low-power and high-power plasma conditions 
are plotted in Fig. 3, for averaged data as previously de- 
fined. The three data sets mostly exhibit the same trends 
and a very good agreement after 1998. Compared with 
the ohmic set, the ECE TS ratio is sometimes slightly 
higher for high-power plasma conditions, reaching 10% 
between 2007 and 2009 and even 15% between 1995 to 
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Pulse No: 35k - 80k. ECE TS ratio 




Timeslices 



FIG. 3. ECE TS ratio against timeslices for the large data 
set, 1995-2009. Ohmic data in black (15k pulses, 135k times- 
lices), low-power in blue ( < Padd < 10 MW, 6k pulses, 
23k timeslices), high-power in red (10 < Padd < 20 MW, 3k 
pulses, 10k timeslices). Years are indicated by dashed lines. 

Pulse No: 70k-80k. ECETS ratio 




Timeslices 



FIG. 4. ECE TS ratio against timeslices for large data set 
limited to C20-C27, 2008-2009, i.e. 5800 pulses, 41k times- 
lices. Ohmic data in black low-power in blue ( 0<Padd < 10 
MW), high-power in red (10 < Padd < 20 MW). Data are 
averaged over 50 timeslices. JET campaigns are indicated. 

1998. A slight offset might be observed between ohmic 
and non-ohmic data. 

Focusing on the ohmic data set, Fig. 2 and Fig. 3 show 
the same variations of the ECE TS ratio, that appear 
to be without clear trends. Large drops of ~20% are 
observed in 1997, 2001, 2003 and 2009. Numerous in- 
creasing or decreasing trends whose amplitudes amount 
to 10-20% can be seen, spanning large numbers of times- 
lices (thousands of them, say typically a few months); 
less sudden variations are also observed. 

Data selection is very broad and spans almost fifteen 
years. ECE TS ratio exhibits many sudden changes; ex- 
plaining the precise causes of most of those changes is 
out of the scope of this paper. For instance, different 
vessel configurations have been used in JET including 
the various divertor, different operation modes, differ- 
ently shaped plasmas etc. It has been checked that these 
changes of divertor configuration do not correspond with 
any obvious evolution of the ratio. 



3. Observations: C20-C27, 2008-2009 

ECE TS ratio shown in Fig. 4 is for a time selection 
limited to the JET campaigns between 2008 and 2009 
(C20-C27, 72000 <= JPN <= 79853). When available, 
HRAD and HRTS data have been selected as well - but 
not shown in Fig. 4. 



Like Fig3, ECE TS ratio shows ±13% variations within 
[0.8-1.05], with the same strong evidence of evolution pat- 
terns well outside a ~8% near-normal scattering. Ohmic, 
low-power and high-power wide sets show almost the 
same evolutions and levels (ohmic and low-power only) 
for the ratio, whereas the high-power ones have a slightly 
higher level in some case, up to 10% for C20. 



4. Key points 

This is one of the main conclusion of the paper: be- 
tween 1995 and 2009, the (temporal) stability of ECE TS 
ratio can not be guaranteed, even for short time-period 
(say a few months), as shown in Fig. 2 and Fig. 3; T e 
measurements of MICH and LIDAR do not compare well, 
even on averaged data. This point is regarded as a prob- 
lematic one from a data coherence point-of-view, as this 
ratio is expected to stay around 1.0 with some reasonable 
margins. 

Large-scale drifts of the ECE TS ratio can obviously 
not been explained satisfactorily by monotonic drifts of 
ECE or LIDAR calibrations. This point will be discussed 
more thoroughly (mostly in (V)). 



D. 1998-2009 reference pulses 

ECE TS ratio variations previously detailed for the 
large data set might be explained by variations of one or 
many plasma parameters. In order to identify possible 
cause- and-effect links, a first step is to select pulses with 
plasma conditions kept as constant as possible. 



1. Description 

As such dedicated reference pulses have been neither 
designed nor performed through the years, a few sets of 
plasma conditions have been chosen and flagged as being 
used as reference pulses: recovery or cleaning pulses with 
a long enough Ohmic phase (3 s). Usually performed for 
the machine to recover from a disruption, they are usually 
made on a daily or weekly basis and do have very similar 
conditions through the years. These recovery pulses were 
not designed to be fully reproducible, e.g. the impurities 
deposits displaced by the disruption they are recovering 
from are expected to affect the plasma composition, but 
no better pulses sets have been found. 

Likewise, references pulses were not designed to main- 
tain a constant T e . 

As detailed in Table II where typical ranges of main 
plasma parameters are listed, 331 reference pulses have 
been manually selected for 58300 < JPN < 79853 (2003- 
2009) corresponding to 5702 timeslices. For dates prior 
to 2003, very few pulses could be found with Bq = 2.41 T 
and I P = 1.96 MA. 
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Some signals are strongly scattered or show non- 
negligible variations: Ohmic power (Pohm), central elec- 
tron density (n e0 ), central electron pressure (P e o), ef- 
fective charge (Zeff), plasma triangularities and T eMAX 
from MICH and LIDAR. Besides T eMAX , the only mea- 
sured signals showing remarkable trends are the vacuum 
field (±0.4%) with similar variations as the total field, 
and Zeff (±25%). 



2. Observations 

As shown in Fig. 5, four diagnostics are available for 
C20-C27 campaigns: MICH, LIDAR, HRAD and HRTS, 
for p = 0.5. Each diagnostic shows similar variations for 
different radial positions p = [0.2,0.3,0.5], even if these 
variations are not always similar from one diagnostic to 
another. Pulse-to-pulse cross-calibration makes HRAD 
follow MICH very closely, as expected. Even if HRTS 
is cross-calibrated against HRAD for the whole C20-C27 
period (see (II A)), variations of both measurements do 
not appear to be linked. Note that LIDAR and HRTS 
exhibit different variations, even if their measurements 
are from the same physics basics. 

As shown in Fig. 5, T e is obviously not constant for the 
reference set: no reference temperature can then be used 
to assess the stability of diagnostics. Similarities between 
the ECE TS ratio and T eMAX MICH and T eMAX lidar 
are to be expected, but depending on the pulse range, 
trends have features sometimes linked to T eMAX MICH 
variations and sometimes T eMAX ,LIDAR : again, no firm 
conclusion can be drawn on the cause of those drifts. 

The ECE TS ratio for reference pulses is plotted in 





Timeslices 

FIG. 5. Electron temperature measurements (keV, p = 0.5) 
against timeslices for reference data set as described in Ta- 
ble II, limited to C20-27 pulses: 124 pulses and 1227 timeslices 
selected. Measurement from different diagnostics are shown: 
MICH, HRAD, LIDAR and HRTS. 



ECE TS ratio 




Signal 






Average 


zb dev. 




in% 


Bo 


[T] 




2.41 


±0.01 




(0.4%) 


Ip 


[MA] 




1.96 


±1.2 x 10" 


-2 


(0.6%) 


Pohm 


[MW] 




0.81 


±0.52 




(65%) 


Wdia 


[MJ] (*) 


0.90 


±5.4 x 10" 


-2 


(6%) 


ne [lxl0 19 m- 3 ] 


2.24 


±0.27 




(12%) 


PeO 


[kPa] 




6.80 


±2.0 




(29%) 


T e M AX,MICH 


[keV] 




2.00 


±0.28 




(14%) 


T e M AX, LIDAR 


[keV] 




2.17 


±0.35 




(16%) 


Zeff 






1.99 


±0.69 




(35%) 


Pn 




(*) 


0.38 


±2.4 x 10" 


-2 


(6%) 


li 




(*) 


1.19 


±3.7 x 10" 


-2 


(3%) 


Q95% 




(*) 


3.67 


±0.21 




(6%) 


plasma volume 


[m 3 ] 


(*) 


80.7 


±1.80 




(2%) 


triangularity (upper) 


(*) 


0.14 


±8.3 x 10" 


-2 


(57%) 


triangularity (lower' 




(*) 


0.18 


±0.07 




(39%) 


elongation 




(*) 


1.56 


±0.12 




(7%) 


midplane centre r-pos [m] 


(*) 


2.94 


±1.1 x 10" 


-2 


(0.3%) 


midplane centre z-pos [m] 


(*) 


0.31 


±2.3 x 10" 


-2 


(7%) 



TABLE II. Plasma parameters for reference data set, limited 
to 58632 < JPN < 79645, 05/03/03 to 14/10/09, 5702 times- 
lices, 331 pulses. 
(★): data computed by EFIT. 



FIG. 6. ECE TS ratio against timeslices for reference data 
set as described in Table II. Years are indicated by dashed 
lines. 

Fig. 6. As in (IIIC2), residues have a near- normal 
(Gaussian-like) distribution similar to the one obtained 
for the large data set - with the mean around 0.95 and 
non- negligible scattering (=half-width) of 0.1. Peak-to- 
peak amplitude is slightly above 30%. 

ECE TS ratio drifts and variations are found to be 
in agreement with the large sets (Fig. 3) and Pearson's 
cross correlation coefficient between ratios from both sets 
is quite high: around 0.7 for JPN > 54500 (from 2003 to 
2009) as well as for C19-C27 (JPN > 70000). 

As stated previously, B and Z E ff are the only mea- 
surements (other than T e ) showing noticeable trends. 
Zeff does not show any variations similar to the ECE 
TS ratio, yet. In Fig. 7, magnetic field (measured by Ro- 
gowskii coils) show large-scale drifts similar to the ones 
observed for ECE TS ratio, but the observed Bo vari- 
ations have a very low amplitude, around ±0.4%; note 
that these experimental variations are below the uncer- 
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FIG. 8. KC1D integrator test against years. Expected values 
= 26214 bits, drift amplitude = 120 bits = ~ 0.2% variations. 



FIG. 7. ECE TS ratio (left y-axis, black) and vacuum mag- 
netic field (Bo, right y-axis, blue) against JPN (top) and 
timeslices (bottom) for reference data set as described in Ta- 
ble II. Years are indicated by lines. 



A possible cause for Bo variations is described and 
analysed in the next section. 



tainty level of KC1D integrator measurement (1.5%, see 
(IID)). 

As stated in (II A) and shown in Fig. 1, MICH and 
LIDAR LOS are not similar, both missing the plasma 
center by ~5 cm according to EFIT reconstruction. For 
reference pulses, r- and z-position are observed to move 
by approximately 2-3 cm at most (see Table II). Note 
that from (IID) these values are believed to be within the 
uncertainties of EFIT reconstruction in the plasma cen- 
ter. Even if those displacements are regarded as real, T e 
profiles in usual JET plasmas as like for the (ohmic) ref- 
erence pulses are not peaked enough to explain this 30% 
peak-to-peak variations. Furthermore, measurements are 
actually volume- averaged for both MICH (usually a few 
cm) and LIDAR (~lcm). 



3. Key points 



IV. OBSERVATION OF SEASONAL PATTERNS 

4. Description 

In JET, magnetic field signals (pick-up coils, loops, 
Rogowskii coils) used to build Bo are all recorded by the 
same type of 8-channel isolated integrator (referred as 
KC1D integrator) used for the principal magnetic dagnos- 
tics. When those KC1D integrators were installed in Au- 
gust 1994 there was some concern as to how the inner 
capacitors would age and perform with respect to tem- 
perature. To estimate the extent of these effects, a very 
stable test-signal generator was permanently installed in 
KC1D and connected to the 8 channels of a KC1D integra- 
tor sample 13 . For each pulse, the data from this module 
are taken and stored as if they were real coil signals. 



For reference pulses, variations of ECE TS ratio appear 
to be mostly disconnected from the variations observed 
on the 18 plasma parameters described in Table II, except 
for the very low amplitude variations of Bo. 

About the similar long-term drifts shown by Bo and 
ECE TS ratio (Fig. 7), variations amplitudes differ from 
more than one order of magnitude: ±0.4% for Bo and 
±15% for ECE TS ratio. Yet long-term trends of both 
quantities show a reasonably good agreement. But Bo 
variations do not affect the computation of T eMAX MICH 
and T e MA x, LIDAR- For ECE diagnostics, frequency- 
position mapping adds radial uncertainties following 
Ar/r -AB/B: ±0.4% for B converts then to ±1-1.5 cm 
at the plasma centre (R~3m), which is below the ex- 
pected uncertainties caused by the midplane mapping. 

From these two points, the expected impact of ±0.4% 
Bo variations on T eMAX are considered as negligible - and 
vice versa. ECE TS ratio and Bq long-term drifts show 
then a good agreement between 2003 and 2009, but no 
direct cause-and-effect link between both can be thought 
of by the authors. 



5. Observations 

The off-plasma signals of KC1D integrator are plotted 
in Fig. 8. Seasonal drifts are observed, as stated in the 
report 13 : "Up until the JET 2000/2001 shutdown, chan- 
nels 2 to 8 show cyclic variation of about ±0.1%, which 
appear to follow winter-summer pattern". The cyclic vari- 
ations observed were considered to be temperature effects 
on these capacitors, but this point has not been studied 
further. 

As plasma measurements of Bq are performed with the 
same type of integrators, they are then expected to suf- 
fer the same seasonal variations. As shown in Fig. 9, 
Bo shows variations quite well matched by KC1D integra- 
tor measurements for the time period mid-2002 to 2009. 
Amplitudes are quite similar: around ±0.4% for Bo and 
around ±0.2% for KC1D integrator drifts. The agreement 
is quite poor for the 1998-2000 time period. 
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Pulses No: 45000-79853. References pulses 
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FIG. 9. Variations of vacuum toroidal field (Bo, left y-axis, 
full symbols) and KC1D integrator drifts (right y-axis, open 
symbols, bits-26000) against JPN (top) and timeslices (bot- 
tom) for references data set as described in Table II. Ad- 
ditional data sets with similar plasma conditions are shown, 
for 45000 < JPN < 58632 with different B values: 2.41 T in 
blue, 2.51 T in magenta, 2.58 T in red. Years are indicated by 
dashed lines. 



Pulse No's 45000-79853. References Pulse's 
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FIG. 11. MICH in-lab measurements (left y-axis, squares) 
and KC1D integrator drifts (right y-axis, bottom, red points) 
against time for 1998-2010. Synthetic signals with a 12-month 
period (black dashed lines) are added on both plots. 



TS ratio above a 1% uncertainty, even for the warmer 
days/months when the night/day temperature variations 
reach up to 15 °C - quite uncommon for Oxfordshire 3 . 
Time constants considered here are consequently more 
likely to be of the order of weeks or months than days. 



V. LONG-TERM STABILITY OF MICHELSON 
INTERFEROMETER 



FIG. 10. ECE TS ratio (left y-axis, black) and KC1D integra- 
tor drifts (right y-axis, blue) against timeslices for references 
data set as described in Fig. 9. KC1D integrator drifts values, 
slightly re-scaled and shifted (offsets < 20 bits), are in red. 
Years are indicated by dashed lines. 



6. Key points 

Given the similarity of Bo and KC1D integrator vari- 
ations, in terms of amplitude and temporal drifts, it is 
likely that Bo variations observed in Fig. 7 and Fig. 9 
after mid-2002 should not been accounted for real events 
taking place in JET plasmas. 

Note that those off- and on-plasma variations affect- 
ing the field measurements are below their experimental 
uncertainty (1.5%, (II D)). 

As stated in (HID) (Fig. 6), the only plasma parame- 
ter showing variations correlated with ECE TS ratio for 
the reference pulses is Bo. The causality of the link be- 
tween (on-plasma) ECE TS ratio and (off-plasma) KC1D 
integrators tests can not be established from the available 
data, or the possibility of a common but still unknown 
cause having similar effects on both signals. Seasonal 
variations of external temperature at JET site are yet 
regarded as the main suspect of the long-term drifts of 
ECE TS ratio after mid-2002 and since 2009. Little can 
be said between 1995 and mid-2002. 

Daily variations of JET site temperature might also 
well trigger similar effects. The corresponding analysis 
has been made and does not show any effect on ECE 



Another way of explaining the observed variations of 
ECE TS ratio is to directly check the variations of T eMICH 
and T eLIDAR , but even for the reference data set that is 
more "constrained" than the large data set, no other T e 
measurement can be used as a reference to assess MICH 
and LIDAR stability. The idea developed in this section is 
to check MICH stability with the in-lab (i.e. off-plasma) 
monitoring performed throughout the years. 



A. In-lab measurements 

In order to check the stability of MICH, in-lab cali- 
brations are routinely performed every few months us- 
ing a calibrated hot source since the early years of JET 
(1984); since 1996, the same hot source (600 °C, built by 
SPECAC in the 90s) has been used. Measurements on 
a cold source (liquid nitrogen) were performed but not 
sufficiently frequently (i.e. monthly). These measure- 
ments were performed by different MICH operators but 
as the set-up for in-lab calibrations is documented and 
reproducible, it is believed that these measurements were 
done in very similar ways. 

Results of MICH in-lab measurements are reported 
in Fig. 11. All the measured interferograms have been 
reprocessed using the same routines and the resulting 



3 http://www.metoffice.gov.uk 
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spectrograms have been integrated over the 80 GHz to 
200 GHz band, which is most frequently used for mea- 
suring T e with the second ECE harmonic in X-mode. 

No cyclic pattern can be easily observed in Fig. 11. 
Data scattering is below ±5% for 1998-2002 and below 
±10% for 2003-2010 which is at least twice lower than 
the scattering of ECE TS ratio. The scattering increase 
of MICH in-lab measurement around 2003 takes place 
approximately at the time where ECE TS ratio starts 
to be well correlated with the seasonal drifts of the off- 
plasma magnetics measurements - as shown in Fig. 10. 
The ECE TS ratio does not show such an increase of 
scattering. 

B. Cryodetector bias voltage 

The bias voltage of MICH's cryodetector is provided 
by batteries delivering a very stable voltage (7 V or 10 V 
depending on the battery model). Sensitivity of the cryo- 
detector to this bias voltage has been characterised. Tak- 
ing 10 V as a reference, detector sensitivity is slightly 
increased (5%) for 7 V and decreased for lower voltages 
(> 10% decrease below 4 V). Sensitivity variations stay 
below 5% for a bias- voltage from 6 to 10 V. This volt- 
age was monitored on a monthly basis and the conse- 
quences of three serious drops (< 6 V) have been care- 
fully examined. The expected impact on T eMAX ,MICH 
is a (possibly slow) > 10% decrease, then a steep recov- 
ery when batteries are recharged. This signature has not 
been observed on JET data; the causality link between 
such uncontrolled variations of the MICH cryodetector 
bias voltage and the drifts of the ECE TS ratio is not 
established then. 



C. Limits 

This monitoring includes the MICH system, the cryo- 
detector and some electronic components. DAQ (hard- 
ware) is different for in-lab and plasma measurements 4 . 
The fifty meters of waveguide, the vessel windows and 
the antenna are not included in this monitoring. 

Components that are included in the monitoring and 
are the more prone to suffer from external temperature 
variations are the cryodetector, which uses liquid helium, 
and the electronic amplifiers, which provided a 50 dB to 
54 dB gain for plasma measurements and 80 dB for in-lab 
measurements and in- vessel calibrations. Their sensitiv- 
ities to temperature variations (~15 °C) have not been 
checked yet. 

Propagation through the waveguide reduces the signal 
by 7 to 15 dB, following an F 2 scaling between 70 GHz 
and 350 GHz. These values are not expected to vary with 



temperature. The in- vessel window certainly slightly 
modifies the EC signal, but should not cause any sea- 
sonal variations. The DAQ's old electronics might be 
sensitive to temperature variations. These three points 
are possibly causing calibration drifts but their ampli- 
tudes are expected to remain quite small, and they are 
not expected to introduce any seasonal variations. 



D. Key points 

This monitoring was historically designed to check the 
time stability of MICH between in-vessel calibrations; 
tracking down such seasonal variations would require at 
least monthly measurements, which is not the case. 

The main conclusion is that the available data from 
MICH monitoring does not show any 12-months cycles. 
Furthermore, data scattering of in-lab MICH calibration 
is a factor of two lower than observed on the ECE TS 
ratio. 

Observed variations of in-lab measurements could 
come from: 

• Experimental errors in the calibration process, like 
hot source position and temperature and different 
transmission in the open waveguide . . . but these 
will not affect T eMICH measurement. 

• MICH calibration drifts that will affect T eMICH mea- 
surement, for example the modification of ageing 
amplifiers ' gain or evolution of the components of 
the cryogenic detector. 

On the other hand, some other effects described in 
(V C) and outside the monitoring scope could modify T e 
measurements. 

Although available data from MICH lab monitoring 
does not support the idea that MICH is mostly respon- 
sible for ECE TS ratio variations, lack of data prevents 
any firm conclusion to be drawn. 

A possibility to definitely track down such calibration 
drifts would be to perform a full in-vessel ECE calibration 
on a monthly basis (at least), during 1 or 2 years. This 
is however not feasible on JET. 



VI. Wth ESTIMATIONS 
A. Description 



A new common DAQ is to be installed for late 2011 or early 2012 



Plasma thermal energy (Wth) is estimated from the 
plasma pressure profiles, i.e. using the kinetic expressions 
following Ref 11 : 
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Pulse No's: 36976-79853. W th ratio with MICH 

Plasma thermal energy for P NB , > 5MW W th ratio with LIDAR Te 




Timeslices 



FIG. 12. Ratio of (kinetic) Wth to (magnetics) Wth, mag with 
MICH (red) and LIDAR (blue) T e profiles in MJ against 
timeslices for references data set. Data selection from JPN 
=36857 to 79853 with P N bi > 5 MW, 672 pulses selected, 2604 
timeslices. Smaller points represent raw data, bigger points 
represent smoothed data. Years are indicated by dashed lines. 

W th =W th , e +W th , i 

W th ,e = \J n e T e dV 
w 7 — Z E ff Ti 

W th ,i ~ =r W th ,e 

6 T e 

Plasma density is taken from LIDAR measurements 
while T e is measured by MICH (W t h,MlCH) and LIDAR 
(W t h,LIDAR)- I° n temperature is measured by charge 
exchange diagnostic 24 . 

Another way to estimate W t h uses magnetics measure- 
ments, defining then W t h,mag- 

W t h,mag WdIA - f Wf ast ,± 

where Wf as t,_L is the stored energy in the perpendicu- 
lar component of the fast particle population, estimated 
from the computation of power deposition of ion cy- 
clotron reasonance heating (ICRH) by PION 8 ' 9 for NBI- 
only plasmas. 

Comparing kinetic W th ,MICH and W t h, LIDAR to the 
diamagnetic W t h,mag provides information on the mea- 
surement accuracy of T e by MICH and LIDAR. Wth 
is then estimated following both ways for purely NBI- 
heated plasmas from the wide set, for Pnbi > 5 MW. 



B. Observations 

Estimations of W th ,MICH/Wth,mag and 
W th ,LIDAR/Wth,mag are shown in Fig. 12. Large 
similarities are observed between both kinetics esti- 
mations and they clearly depart from the diamagnetic 
one, with variations between -40% and 10%. Most of 
the time W th ,MICH < W t h, LIDAR, coherent with the 
observations of ECE TS ratio shown in Fig. 2 and Fig. 3: 
consequently Wth, LIDAR estimation is generally slightly 
closer to the diamagnetic Wth than W t h,MlCH- No 
other systematic trend is observed. 



The same exercise has been done for three other cases: 
first with purely Ohmic plasma where Wth = Wdia on 
a larger pulses database, and then with a condition on 
plasma central density for the NBI-only plasmas (central 
n e0 < 5 x 10 19 m~ 3 and n e0 > 5 x 10 19 m~ 3 ) to check 
the impact of density on the Ti measurements. These 
additional cases lead to very similar conclusions as the 
ones drawn for the NBI-only case. 



C. Key points 

The comparison of kinetic or diamagnetic estimations 
of W t h does not provide any definitive information on 
the accuracy of MICH or LIDAR measurements. It seems 
that MICH is slightly under-estimating T e , but both of 
them lead to a W t h generally 20-30% below the diamag- 
netic Wth, for both Ohmic and NBI-only plasmas. 

VII. ECE TS RADIAL DISCREPANCIES 
A. Description 

Relatively low signal-to-noise ratio and lower radial 
definition of LIDAR T e profiles make it impossible to 
satisfyingly track down ECE TS position errors. Radi- 
ally well-defined profiles measured by HRTS are routinely 
available since April 2008 (JPN -72000, from C20). Af- 
ter midplane mapping, they can be compared to the high 
time and high (radial) resolution HRAD profiles. 

A long lasting issue within the JET community is that 
HRAD and HRTS profiles are not always similar (e.g. 
as described in Ref 2 ' 16 ); these mismatches are referred 
to as ECE TS radial discrepancies. As the radial posi- 
tion of HRTS profiles comes from an independent calibra- 
tion while ECE profiles requires a radial-frequency map- 
ping relying on magnetic reconstruction (as described in 
(II D)), radial corrections are usually applied on the ECE 
side. Profiles features such as the T e pedestal in H-mode 
indicates that at least part of these mismatches can often 
be corrected by a radial shift or a slight modification of 
the vacuum magnetic field (< 1%, less thant the mea- 
surement uncertainty), following ~/c corr B and then ap- 
plying a non-constant radial shift. Most of these points 
are carefully addressed in Ref 16 . 

Observables for ECE TS radial discrepancies are de- 
signed in order to estimate the radial differences between 
HRAD and HRTS profiles. Fig. 13 shows how this esti- 
mation is done for r=3.5m for a typical JET pulse. The 
situation is the same with p profiles. Plasma centre (p~0) 
is usually missed by both diagnostics and the edge brings 
no additional information in ohmic and L-mode; the fo- 
cus is then arbitrarily set to the gradient zone following 
median(p H RAD - Phrts,0.3 < p < 0.8). 

The correction of the ECE TS radial discrepancies, and 
the various ways to do so, have been intensively discussed 
in the JET community for years. No definite answer has 
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Pulse No: 75000 - 1 = (9.5s-1 0.5s) 
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FIG. 13. ECE TS radial discrepancies for typical JET plasma 
(JPN 75000, t=9.5-10.5s) for r=3.5 m. Plasma centre and sep- 
aratrix are respectively at ~3m and ~3.9m. Left: HRAD 
[red] and HRTS [blue] T e profiles; data (points), average 
(lines). Centre: Deviations from time-average for each di- 
agnostic. Right: estimation of Thrad — t H rts for r=3.5m 
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FIG. 14. ECE TS ratio (left y-axis, black points) and ECE TS 
radial discrepancies median(p H RAD — Phrts,0.3 < p < 0.8) 
(right y-axis, p units in %) against timeslices for references 
pulses for C20-C27, where HRAD and HRTS data are avail- 
able. Ohmic conditions, 147 pulses, 2878 timeslices. ECE 
TS radial discrepancies are non-corrected (cross, blue) and 
corrected (square, red) of the ECE TS ratio variations. 
Month/years are indicated by dashed lines. 



been found yet, because such corrections must ensure 
that HRAD and HRTS profiles match at the centre, in 
the gradient zone and at the plasma edge. Note that 
the situation worsens because the degree of the ECE TS 
mismatch is not always the same for ohmic, L- and H- 
mode plasmas. 

In this part the possibility of a link between the ECE 
TS ratio (from MICH and LIDAR) and the ECE TS radial 
discrepancies (from HRTS and HRAD) is examined. 



B. Observations 



In Fig. 14, profiles mismatches are expressed as am- 
plitude mismatches (ECE TS ratio, made of T eMAX mea- 
surements from MICH and LIDAR) and radial ones (ECE 
TS radial discrepancies, averaged over 0.3 < p < 0.8) for 
C20-C27 reference pulses where data from all these di- 
agnostics is available. Only ohmic plasmas are analysed. 
Both quantities are well correlated. Radial discrepancies 
are mostly ranging from —0.1 to (p units). 

Correcting HRAD T e profiles by increasing their values 




p units (%) 



FIG. 15. Distribution of ECE TS radial discrepancies (p units 
in %), non-corrected (blue) and corrected (red) of the ECE TS 
ratio variations. Same pulses/timeslices selection as Fig. 14. 



according to the ECE TS ratio, i.e. multiplying T eHRAD 
by T eMAx L iDAR/ T e M Ax,MICH, leads to a decrease of 
the ECE TS radial discrepancies to a -5% to 5% range 
(Fig. 14, with Fig. 15 for the distribution plot). This cor- 
rection acts like a re-centring of the radial discrepancies 
around 0%; data dispersion is not reduced yet. 



C. Key points 

ECE TS radial discrepancies might be explained as an 
amplitude shift following the ECE TS ratio - at least 
partially. Taking into account the variations of the ECE 
TS ratio is regarded as a first-step correction. 

Cross-calibrations of HRTS and HRAD against MICH 
emphasises the role of MICH, as any possible drifts of 
MICH calibrations will appear on both ECE TS ratio 
and radial discreprancies sides. Possible drifts of LIDAR 
calibration appear only on the ratio side. Differences be- 
tween corrected radial discrepancies and ratio variations 
are then mostly showing different behaviours between LI- 
DAR central T e and HRTS T e in the gradient zone. 



VIII. DISCUSSION & CONSEQUENCES 

A. Long-term stability of JET T e measurement 

As defined in (II), the T eMAX quantity is a useful means 
of estimating the central electron temperature of JET 
plasmas and of comparing T e measurements done by di- 
agnostics with different lines of sight and/or diagnostics 
that do not always access the plasma centre. 

This paper's main finding is that the observed 15 
years of central T e measurements (1995-2009) featured 
in Fig. 3 and Fig. 4 mostly range between 0.5 and 1.25 
(mean around 0.95) and exhibit non negligible variations 
through the years, with a 30% uncertainty observed on 
non-averaged ratio (see (IIIC)). Large-scale structures 
and patterns are identified as large-scale temporal vari- 
ations, amounting up to ±20%. The complexity of the 
observed time drifts does not call for a simple explana- 
tion. Such variations are observed as well on more recent 
JET campaigns (C20-C27) in 2008 and 2009 (IIIC 3). 
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No measurement for T e could be used to assess the 
time stability of any of the studied diagnostics (IIIC, 
HID), even for the C20-C27 JET campaigns when four 
diagnostics were available - only three being calibrated in 
a way that their relative time variations are actually inde- 
pendent: MICH, LIDAR and HRTS. Analysis is therefore 
limited to relative variations between different measure- 
ments. 

Following (HID), variations of ECE TS ratio are dis- 
connected from the 18 observed plasma parameters listed 
in Table II for carefuly selected reference pulses, apart 
from very low variations (±0.4%) of the magnetic field 
B . 

They appear yet to be quite satisfyingly correlated 
with the 12-month-period patterns shown by off-plasma 
KC1D integrator signals for the 2002-2009 time period (IV, 
Fig. 10). Those variations of KC1D integrator are believed 
to be caused by external or ambient temperature. 

The analysis of the in-lab monitoring of MICH (VA), 
consists of measurements against the same calibrated mi- 
crowave source performed every few months since 1996, 
in order to check the long-term stability of the diagnostic. 
It is not possible to draw any firm conclusion on the link 
between drifts of ECE TS ratio and in-lab measurements 
performed on MICH, except that this monitoring shows 
that drifts of the ECE TS ratio cannot be explained only 
by the MICH drifts observed in-lab. 

Regular calibrations of LIDAR are regularly per- 
formed: measurements of spectral functions and detector 
quantum efficiency are required to calculate plasma's T e . 
Such calibrations have been done at random times of the 
year, since at least 15 years. These calibrations have al- 
ways been very stable, so the authors would not expect 
any temperature dependance in them. 

Off-plasma monitorings of both MICH and LIDAR do 
not show therefore any drifts that have been linked with 
the observed seasonal drifts of the ECE TS ratio. 

For the 1995-2009 time period, JET T e measured by 
ECE and TS diagnostics show large-scale slow-varying 
time-dependent variations up to ±20%; calibrations 
drifts are suspected, and quite clear correlations are 
found with seasonal variations of the JET ambient tem- 
perature since mid-2002 and 2009. 

The authors are not aware of such seasonal drifts be- 
tween the ECE and TS systems installed and operated in 
any other large sized tokamak. As JET systems appear 
to suffer such sensitivity at an unexpected scale, it would 
be recommended to perform similar systematic checks on 
the electron temperature measurements on other large- 
sized tokamaks, if absolutely calibrated ECE and TS sys- 
tems are available. 



B. Revised uncertainty for overall measurement of JET 
electron temperature 

As long as the cause (s) of the observed variations of 
ECE TS ratio remain (s) unknown, it should be consid- 



ered as uncontrolled uncertainties affecting JET T e mea- 
surements. The next paragraphs propose a revised un- 
certainty of JET T e measurements from the experimental 
observations presented in this paper, i.e. the 1995-2009 
time period. Future measurements are very likely to be 
similarly affected as well. It focuses on central T e ; no 
better uncertainty is expected for non-central measure- 
ments. 

Experimental uncertainty can be estimated to re- 
sult from a systematic, non-time dependent uncertainty 
(5T/T) ca i accounting for absolute calibration and a time- 
dependent uncertainty (ST/T) t caused by drifts, varia- 
tions or just plain noise, following: 

ST/T ~ yVr/T)cai 2 + (<JT/T) t 2 

Estimated uncertainties of MICH and LIDAR calibra- 
tions i.e. (5T/T) ca i are within 5-10% (II A). Experimen- 
tal large-scale observations reported in this paper repre- 
sent (5T/T) t and are estimated at 15-20%, leading to the 
following estimation for the overall uncertainty for JET 
central T e measurement: 

ST/T - 16 - 22% 

For non-central regions, issues specific to ECE ra- 
dial localisation arise, reducing therefore the accuracy 
of MICH measurements. LIDAR measurements are not 
affected. 

HRAD is cross-calibrated against MICH for each pulse, 
so T eHRAD uncertainty can not be lower than MICH data 
in the plasma center. Same remark for non-central T e . 

Absolute calibration of HRTS relies on HRAD for C20- 
C27, therefore, for this pulse range, it also relies on 
MICH. As long as no MICH-independent calibration is 
used, uncertainties of central and non-central T eHRTS can 
not be lower then MICH uncertainties. 

Possible drifts of HRTS calibration are not assessed in 
this paper, but it can be observed in Fig. 5 that variations 
of HRTS measurements for p = 0.5 might be different to 
those of MICH or LIDAR. Conclusions for p = 0.2 and 
0.3 are similar. The uncertainty of HRTS data can not be 
estimated to be lower than the estimation for the three 
other diagnostics. 

It is worth stressing that the confidence intervals re- 
sulting from absolute calibrations of MICH (±10%) and 
LIDAR (±5%) - see (II A) - are clearly dwarfed by 
the time-dependent drifts observed, which are possibly 
caused by drifts of MICH and LIDAR calibrations. 



C. Future work 

Using MICH, LIDAR and HRTS as absolutely cali- 
brated diagnostics, providing independent measurement 
of the same quantity (T e ), it is hoped that enough data 
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could be gathered, since C20 (April 2008), in order to dis- 
criminate which diagnostics suffer from calibration drifts. 
Few things can be done for the time period before 2008, 
where MICH and LIDAR were the only calibrated diag- 
nostics. 

For the next years of JET operation (2011 and after), 
the authors recommend a reference pulse to be designed 
and performed on a weekly basis at least for, say, 2 to 3 
years. It is of course not possible to guarantee the per- 
fect stability of the plasma parameters, but methodically 
repeating the same pulse will certainly help a lot to track 
the stability of the T e measurements. 

In parallel, in-lab monitoring of the diagnostics should 
be performed when possible. The frequency of the in- 
lab monitoring of MICH will be increased to a weekly or 
twice per month period. Ways of monitoring LIDAR and 
HRTS systems should be investigated. 

IX. CONCLUSION 

Between 1995 and 2009, measurement of JET T e has 
been carefuly reviewed using the two main diagnostics 
available for this time period: MICH (ECE), LIDAR 
(TS). The overall uncertainty of T e measurements, based 
on these two diagnostics, is estimated to be as high as 16- 
22%, mostly made of large-scale time drifts. 

From mid-2002 to 2009, the ratio of temperature mea- 
sured by the two diagnostics (namely the ECE TS ratio) 
shows clear patterns with a 12-months period and a 15- 
20% amplitude. These drifts are regarded as calibration 
drifts, possibly caused by unexpected sensitivity to un- 
known parameters. External temperature on the JET 
site might is the best parameter suspected so far. The 
final cause of such variations has not been found, nor 
the way it affects MICH, LIDAR and/or (possibly) HRTS 
measurements. 

It is not clear whether it is direct impact on one of 
the diagnostics or indirect via magnetics measurements, 
which have been shown to have seasonal drifts. 

Solutions aiming at tracking down these unexpected 
uncertainties of JET T e have been detailed and can be 
performed during the following JET campaigns (C28+, 
after October 2011), including reference pulses designed 
to be highly reproducible over time (new ITER-like wall 
permitting) and specific off-plasma monitoring of the di- 
agnostics - when applicable. 

Whatever causes these JET T e drifts, this experimen- 
tal issue is regarded as crucial for data quality. Long-term 
observed ECE TS radial discrepancies have been shown 
to be linked to these seasonal drifts. Tackling this issue 
and elucidating the cause of these unexplained seasonal 
variations is of great importance, impacting to the qual- 
ity of JET T e measurements. 
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